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ABSTRACT 

We have used a set of dedicated astrometric data from the Hubble Space Telescope to measure the 
parallax and proper motion of the nearby neutron star RX J0720.4— 3125. At each of eight epochs 
over two years, we used the High Resolution Camera of the Advanced Camera for Surveys to measure 
the position of the B = 26.6 target to a precision of ~ 2 mas (~ 0.07 pix) relative to 22 other stars. 
From these data we measure a parallax of 7r = 2.8 ± 0.9 mas (for a distance of 360^qq pc) and a 
proper motion of fi — 107.8 ± 1.2 mas yr _1 . Exhaustive testing of every stage of our analysis suggests 
that it is robust, with a maximum systematic uncertainty on the parallax of 0.4 mas. The distance 
is compatible with earlier estimates made from scaling the optical emission of RX J0720.4— 3125 
relative to the even closer neutron star RX J1856.5— 3754. The distance and proper motion imply a 
transverse velocity of 180^40 km s_1 > comparable to velocities observed for radio pulsars. The speed 
and direction suggest an origin for RX J0720.4— 3125 in the Trumpler 10 OB association ~ 0.7 Myr 
ago, with a possible range of 0.5-1.0 Myr given by the uncertainty in the distance. 
Subject headings: astrometry — stars: individual: alphanumeric: RX J0720.4— 3125 — stars: neutron 
- X-rays: individual (RX J0720.4-3125) 



1. INTRODUCTION 

One of th e many interesting results from ROSAT All- 
Sky Survey (|Voges et al.lll999f ) was the discovery of seven 
objects that appear to be nearby, thermally-emitting 
neutron stars that have little if any magnetq s pheric 
emission (for recent reviews , see lHaberll 120041 l200l 
Ivan Kerkwijk &: Kaplan! l2006f ) . These objects, known 
most commonly as "isolated neutron stars," (INS) are 
distinguished by their long spin periods (> 3 s, when 
measured), largely thermal spectra with cool tempera- 
tures (kT < 100 eV), faint optical counterparts, and lack 
of radio emission. 

Thermally-emitting neutron stars have been the tar- 
gets of many observations, as they can potentially be 
used to constrain the equation of state (EOS) of neu- 
tron stars, and thereby explore nucl ear physics in realms 
inacc essible from laboratories (e.g.. lLattimer fc Prakasr] 
2000). Two main approaches are used. The first, us- 
ing the spectrum, seems simple: determine the effective 
angular size from spectral fits, multiply by the distance 
(obtained by other means), and one has the apparent 
radius. This radius can be converted into the physi- 
cal radius through use of mass. The radius is the cru- 
cial quantity in differentiating between EOS, as most 
EOS predict a distinctive but small range of radii for 
a large range of masses. The second method is to use 
measurements of an ensemble of neutron stars to con- 
strain cooling curves, which are themselves sensitive to 
the EOS and the interi or composition (for a review, 
lYakovlev fc Pe thick 2004) . This also seems simple: mea- 
sure the received flux of an object, convert to a lumi- 
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TABLE 1 

HST ACS/HRC F475W Observation Summary 





Date 




Orientation 51 


TO475 b 


Epoch 


CUT) 


(MJD) 


(dog) 


(mag) 


1 


2002 Jul 04 


52460.0 


353.0108 


26.66 


2 


2002 Sep 15 


52532.6 


72.5457 


26.54 


3 


2003 Jan 06 


52645.1 


175.9862 


26.45 


1 


2003 Mar 16 


52714.0 


254.8005 


26.57 


5 


2003 Aug 04 


52855.1 


28.0808 


26.39 


6 


2003 Sep 21 


52903.0 


77.1792 


26.45 


7 


2004 Jan 07 


53011.9 


178.4808 


26.42 


8 


2004 Mar 19 


53083.1 


258.1951 


26.64 



Note. — For each epoch, eight exposures were taken 
during two orbits. For epochs 1-4, the integration times 
for exposures 1—4 were 600 s, while those for exposures 5—8 
were 630 s. For epochs 5—8, the exposure times were 625 s 
in the first orbit, and 660 s in the second. The dither offsets 
along the x and y axes of the detector (written in base eight, 
where, e.g., 2.2 is 2+§), were [0.0, 0.0], [2.2, 2.2], [4.4, 4.4], 
[6.6, 6.6], [1.1, 5.5], [3.3, 7.7], [5.5, 1.1], and [7.7, 3.3]). 
a This is the value of the orientat keyword, roughly cor- 
responding to the j/-axis on the HRC. The transformation 
angles in Table [5] are relative to these. 

b The F475W magnitude of RX J0720.4-3125, measured 
at each epoch. The uncertainty on each measurement is 
~0.09 mag. 



nosity using the distance, and with the age place that 
object on a cooling diagram. The closest of the INS, 
RX J1856.5— 3754, has been t he subject of much inquiry 
for just these purposes (e.g.. iWalter fc Lattimerl 120021 : 
iDrake et al1l2002t iBraie fc Romanill2002f) . 

While the basic approaches — inferring radii from 
broad-band spectral fits and comparing cooling lumi- 
nosities to models — are easily stated, complications 
arise from uncertainties in distances and ages; addi- 
tional complications may result from the presence of 
strong surface magnetic fields and unknown surface 
compositions. Progress on these last issues has been 
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Fig. 1. — Combined image of all 64 HRC exposures of the field of RX J0720.4— 3125 (made using multidrizzle; IKoekemoer et alj|2002ft . 
The objects used for the astrometry are labeled (remaining objects were either too faint or too extended to be used). We also show the 
fitted proper motions of those objects, with the ellipses indicating the 1<t uncertainties on the proper motion. Note that the proper motion 
of RX J0720.4— 3125 (object "X") has been scaled down by a factor of 1/10 compared to the other objects in the field. The scales and 
orientation are indicated in the upper left corner. 



made recently from X-ray spec tr oscopy (lHaberl et al.l 
20031: Ivan Kerkwiik et all [2004; H aberl et alTl2 004b; 



Zane et al.1 120051 ) and timing ( Kaplan &; van Kerkwijkl 



2005a( bT). but distances and ages are difficult to mea- 
sure. The technique necessary for both is high- 
precision optical astrometry, as shown by the example of 
RX J1856.5— 3754: with a series of Hubble Space Tele- 
scope (HST) observations it was possible to measure 
the parall ax dKaplan. van Kerkwiik. &; Ander son 2002b, 
hereafter |KvKA02t iWalter & Lattimerl I2002D . and the 
proper motion traced th e object to a n OB association 
where it was likely born (|Walter|[2001l ). 

While RX J1856. 5-3754 is the brightest and closest 
of the INS, and has thus garnered most of the atten- 
tion, it is best to try to use multiple sources. This es- 
pecially since each source appears to have peculiarities, 
be it stronger or weaker timing noise, stronger or weaker 
features in the X-ray spectra, presence or absence of long- 



term variations, or the prese nce of an Ha nebulae (in 
the c ase of RX J1856. 5-3754: Ivan Kerkwiik fc Kulkarnil 
[200lh . Likely, secure results will only be obtained if we 
understand and can correct for these differences. Here, 
we discuss high-precision astrometric observations of the 

second brightest source, RX J0720.4— 3125. 

RX J0720.4-3125 was discovered by lHaberl et all 
(1997) as a soft (kT ~ 80 eV), bright X-ray source 
in the ROSAT All-Sky Survey. Given its low hydro- 
gen column density (Njj ~ 1 x 10 20 cm -2 ), nearly 
sinusoidal 8.39-s pulsations, relatively constant X-ray 
flux, and faint [B = 26.6 mag ) , blu e optical counter- 
part (jKulkarni fc van Kerkwiik! [l99l iMotch fc Haberll 
1998), it was classified as a nearby, isolated, thermally- 
emitting neutron star. X-ra y timing observations 
(|Kaplan fc van Kerkwiik] l2005al) give a characteristic 
age of 2 Myr and a magnetic field strength of 2.4 x 
10 13 G, consistent with suggestions of the source be- 
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ing an off-beam, moderately strong-fiel d radio pulsar 
(jZane et al.l l2002t Kaplan et all l2002af ). While the 
above properties make RX J0720.4— 3125 a prototyp- 
ical INS, X-ray monitoring over the last few years 
has uncovered unique behavior: the X -ray spectrum 
and pulse shape have been evolving (|de Vries et alJ 
[200l IVinket al.ll200l. This may indicate free pre- 
cession ( Haberl et al.ll2006f). although this inte rpretation 
is not unique ( jvan Kerkwiik fc Kaplan! 12006). Scaling 
to RX J18 56.5— 3754 , its distance was est imated to be 
~ 300 pc (|KvKA02fl . iMotch et all (|2005l ) measured a 
proper motion of 97 ± 12 mas yr _1 from ground-based 
astromctry, and suggested the source might originate in 
the OB association Trumpler 10, something we inde- 
pendently der ived from part of the data presented here 
(|Kaplanll200l . 

In this paper, we present a detailed analysis of HST 
observations of RX J0720.4— 3125 specifically designed 
for astrometry, which we use to measure its proper mo- 
tion and parallax. The organization of this paper is as 
follows. First, in § [SJ we present the astrometric obser- 
vations and discuss additional observations we used for 
photometry and for tying our astrometry to the Inter- 
national Coordinate Reference System (ICRS). We also 
discuss the basic reduction of these data. Then, in § [3l 
we use the photometry to determine photometric paral- 
laxes of the reference stars, which we use in § |4] for the 
parallax. We discuss the implications in §[5]and conclude 
in § H 

Our analysis is complex and necessarily involves 
choices among alternate schemes. Previous attempts 
to measure parallaxes with H ST made choices that 
were later called into qu estion (jWalterl 12001 lKvKA02l ; 
I Walter fc Lattimerl l2002j ). To clarify our analysis, and 
to set out the framework for future work, we give full de- 
tails on the parallax measurement in App. [X] We then 
show that our results are robust to the choices we made 
by describing various alternate analyses in App. [Bj In 
what follows, we define our proper motions in Right As- 
cension and Declination (/i Q , /is) such that the scales are 
the same and no cos<5 term is necessary. All uncertain- 
ties are la unless otherwise indicated. In addition to the 
results presented in Tables 0H6l we make the raw data 
available electronically. 

2. OBSERVATIONS, REDUCTION, AND ANALYSIS 

Our main data set consist of 64 exposures taken dur- 
ing eight visits of the field of RX J0720.4-3125 with 
HST; a summary is given in Table [TJ All data were 
taken with the High Resolution Camera (HRC, with a 
plate-scale of 28.27 mas pix -1 ) and the F475W (SDSS 
g') filter, which has good sensitivity to the blue colors of 
RX J0720.4-3125. 

We tried to optimize our observations for measuring an 
accurate parallax in three ways. First, we observed over 
two years, to minimize the covariance between proper 
motion and parallax, and verify repeatability. Second, 
we observed four times per year rather than just at the 
parallactic maxima, to reduce possible systematic effects 
of observing at orientations different by 180° only. Third, 
for the eight exposures in each epoch, we chose dither 
positions optimized for sampling fractional pixel phase. 

In our analysis, we also use three sets of observations 
presented previously. For photometry, these are B- and 
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111475-R, 111475-11150+1.5 

Fig. 2. — Color-color relations used to infer B — R and R for stars 
for which the ground-based photometry was inaccurate. Shown 
are relations between B — R and 771475 — R (circles), B — R and 
TI475 — TI50 (squares), and 0150 — -Rand 01475—0150 (lozenges), with 
the 01475 — m,5o magnitudes offset by 1.5 mag for clarity; the zero 
points are indicated vertical dotted lines. Overdrawn (long-dashed 
lines) are the empirical quadratic fits from Table [3] to all points 
except RX J0720.4— 3125 (crossed points), as well as the relation 
between B — R and 01475 — R ex pected from synthetic photometry 
(dotted line: ISirianni et al.ll2005fi . 

R-band images taken with th e Low Resolution Imaging 
Spect rometer (LRIS) on Keck (Kulka rni fc van Kerkwiikl 
1998), and images taken with the Sp ace Telescope 
Imag ing Spectrograph (STIS) on HST ([Kaplan et al.l 
2003), and for our astrometric tie, B-band observa- 
tions taking with the Focal Reducer, low-dispersion 
Spectrograph (FORS1) on th e Very Large Telescope 
(jMotch. Zavlin. fc HabeTl 12001 ). 

Below, we describe first the reduction of all data sets, 
and then the way we measured positions from the HRC 
images, tied our astrometry to the International Celestial 
Reference System, and obtained photometry from the 
HRC, STIS and LRIS data. For reference, we show a 
combined image of all HRC data in Fig. [TJ with stars 
used in our analysis labeled. 

2.1. Reduction 

For each data set, the basic reduction was standard: 
bias and dark current were subtracted, and possible 
pixel-to-pixel sensitivity variations were corrected for us- 
ing flat fields. For the HST data, these steps are taken 
in the standard pi pelines, while for the LRIS data they 
were performed bv lKulkarni fc van Kerkwiildlfl99l . For 
the FORS images, which we retrieved from the archive 
and re-analyzed, we determined the bias level from the 
overscan regions and constructed a flat field from dawn 
sky images. 

In preparation for further analysis, we flagged bad 
pixels and made averages. For the HRC data, we 
used a procedure sim ilar to that of multidrizzle 
(|Koekemoer et all 12002). but which does not correct for 
distortion. First, we removed pixels flagged as bad in the 
data quality array, as well as two obvious bad columns 
[upwards from pixel (x,y) = (200,548) and (617,228)]. 
Next, we subtracted the sky level and resampled on an 
eight times finer grid, for which the dither positions cor- 
respond to integer pixel offsets. Using these offsets, we 
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TABLE 2 

Photometry and photometric parallaxes for stars in the field of RX J0720.4— 3125 







Observed 








. Inferred . 










B 


"1475 




R 




(B - fl) ad 






""phot 


ID 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


Source 


(mag) 


(mas) 


X 


26 62H 71 


26.51(4) 


26.45(5) 


27.0(4) 


27.0(4) 


_ n 4(41 


B - 


R 






A 


20 51 7(41 


20.221(5) 


20.150(10) 


1 9 41 8(31 


1Q 418(31 


1 099(^1 


B - 


R 


•1.5 


0.105 


rj 


20 789(41 


20.468(5) 


20.337(10) 


19 61 2(31 


1 9 61 2(31 


1.177(5) 


B - 


R 


•1.8 


0.110 


104 


23 60(21 


23.012(5) 


22.272(11) 


21 425(1 5"! 


21 425(1 51 


2.17(3) 


B - 


R 


7.5 


0.166 


105 


27.0(2) 


25.80l(l8) 


24.578(15) 


23 84(21 


23 84(21 


2 68( 51 




- R 


9.5 


0.134 


106 


26 47(1 3"! 


25.653(18) 


z4.o91(14) 


23 66(21 


23 66(21 


2.73(5) 


772475 


- R 


9.7 


0.165 


107 


23 384(91 


22.920(6) 


22.448(11) 


21 600(41 


21 600(41 

ii _1_ . UVIU1 T ) 


1 784(101 


B - 


R 


6.6 


0.101 


108 




26.77(4) 


25.50(2) 


24 98(81 


24 82(91 


2.67(8) 


777475 


■ "150 


9.4 


0.082 


1 no 




24.410(12) 


23.499(12) 






9 A4?3\ 


1 1 1 ;■, ' -) 


- R 


0.0 


n 1 An 


110 




22.094(3) 


21.022(10) 


90 1 


90 1 Q^?^ 




B - 


R 


8.7 


51 


111 


23 731 f 1 9 s ) 


23.235(8) 


22.698(11) 


21 £13f4^ 


91 813?4^ 


1 Q1 8?1 2^ 

_L . i7 _1_ CM J- 1 


B - 


R 


6.9 


0.105 


112 




26.43(3) 


25.127(20) 


24.4(2) 


24.47(7) 


2.69(7) 


777475 - 


■ m 50 


9.5 


0.102 


113 




26.79(4) 


25.60(3) 




24.86(9) 


2.62(9) 


777475 - 


- "150 


9.1 


0.072 


111 


22 514(61 


22.024(5) 


21.404(10) 


20 514(31 


20 514(3*1 


2 000(61 


B - 


R 


7.1 


0.21 


115 


21.774(4) 


21.392(6) 


21.190(10) 


20.377(3) 


20.377(3) 


1.398(5) 


B - 


R 


5.6 


0.111 


116 


23.045(8) 


22.422(8) 


21.441(10) 


20.569(4) 


20.569(4) 


2.476(9) 


B - 


R 


8.5 


0.38 


117 


27.1(3) 


26.43(3) 


25.30(2) 


24.48(6) 


24.53(7) 


2.59(7) 


777475 - 


■ "150 


9.0 


0.077 


118 


20.510(4) 


20.225(5) 


20.179(10) 


19.498(3) 


19.498(3) 


1.012(5) 


B - 


R 


1.2 


0.087 


119 


19.689(3) 


19.461(7) 




18.855(3) 


18.855(3) 


0.834(4) 


B - 


R 


3.5 


0.087 


120 


24.73(3) 


23.965(11) 


22.507(11) 


21.997(8) 


21.997(8) 


2.73(3) 


B - 


R 


9.8 


0.36 


121 


25.86(7) 


25.115(17) 


24.036(13) 


23.193(15) 


23.193(15) 


2.62(4) 


"M75 


- R 


9.1 


0.154 


122 


24.135(16) 


23.459(5) 


22.151(11) 


21.519(6) 


21.519(6) 


2.616(17) 


B - 


R 


9.1 


0.33 


123 


25.92(16) 


25.229(14) 


23.861(12) 


23.22(3) 


23.25(5) 


2.71(5) 


"1475 - 


■ "150 


9.7 


0.191 


121 


27.3(3) 


26.304(20) 


24.759(16) 


24.56(7) 


24.31(6) 


2.75(6) 


777475 - 


- 7"50 


9.9 


0.132 



Note. — Observed: All errors quoted exclude the zero-point uncertainties to which the me asur ements are on 
the Vega system. These are < 0.02 mag for all but "150, which is only roughly calibrated (see c12.4|l . Reasons for 
missing measurements are: 108 and 112 (B): too faint; 113 (B and R): too close to star A; 119 ("750): overexposed. 
Inferred: ij a( j and (B — -R) ac j are our adapted values, inferred from the magnitudes listed under Source. The inferred 
absolute magnitudes Mu ms have uncertainties around 0.4 mag, and the corresponding uncertainties in the photometric 
parallaxes n p ^ ot are about 20% (see J3]l- 



TABLE 3 

Color transformations used to infer B — R and R for stars in the field. 



Color 




Empirical Relation 




rms 
(mag) 


B - R 
B - R 

777 5 - R 


= +0.116 
= +0.915 
= +0.654 


+1.173 (ni475 - R) +0.067 ("1475 
+2.270 (m 4 75 - rn 50 ) -0.700 ("1475 
+0.759 ("1475 - "150 ) -0.579 ("1475 


-R) 2 

- "»50) 2 

- "150) 2 


+0.030 
+0.049 
+0.028 



registered the images for each epoch and constructed an 
initial guess at a cosmic-ray free, average image from the 
two lowest values at each position. 

For each image, we flagged as cosmic-ray hits all pixels 
in excess by more than 5a of those in the initial aver- 
age (resampled back to the original resolution) as well 
as all adjacent pixels. Here, we estimated a by adding, 
in quadrature, the read-out noise, the Poisson noise due 
to sky and signal (as inferred from the average), as well 
as terms equal to 5% of the signal and 5% of the range 
in signal in the surrounding 3x3 pixel box. The latter 
terms ensure centers of bright stars are not incorrectly 
flagged as bad due to registration errors, focus changes, 
etc. We also constructed a final average for each epoch 
from those cleaned images. 

For the STIS data, the analysis was similar, except that 
no resampling was necessary, since integer pixel offsets 
were used. For the LRIS and FORS images, we used 
more standard methods to correct for cosmic-ray hits, in 
which they are identified by their narrow width compared 



to the seeing, and replaced by interpolation. For the 
photometry and astrometry below, averages were made 
of these filtered images, registered to integer pixel shifts. 

2.2. Position Measurements 

We determined positions using the "effective" point- 
spread function (ePSF) technique of lAnderson fc Kind 
(2000). We use the ePSF determined for the HRC with 
the F4 75W filter by lAnderson k Kind (|2004L hereafter 
lAKOl . and use a variation of their procedure for the fits. 
In App. IB. 11 we describe where we deviate and discuss 
further variations; fortunately, our results do not depend 
much on our choices, although the method described here 
gives slightly superior precision. 

In our procedure, we fit stellar images on each expo- 
sure using the pixels in a radius of 2.56 pixels around an 
initial guess for the centroid (determined from a Gaus- 
sian fit in the average image), weighting data points 
with uncertainties a = (RN 2 + max(.F, sky)/G) 1 / 2 (where 
RN — 6 DN is the read-noise, which we determined 



The Distance of RX J0720.4-3125 



5 



22 

ttj 



24 



26 









Vl 1 | 1 1 -J, 1 | 1 1 1 1 | 1 


- 




*\ 


\ m V \ 

S \ 0110 

\ \ u O \ 
\^ \^ IW ^ \ 


- 






\ 13111 w 02 \ 

\ \ e\ 
s \ 120 * 
\ \ e \ 

\ X \Cl©@-i23 

\ ^ ^ 








\ - 

\ ra \ 




"+7 


i i 


"ill 



J i i i i I i i i i I i i i i L 

12 3 

B-R 



Fig. 3. — Color-magnitude diagram, showing R magnitude versus 
B—R color for the objects in Table|2] We also show main sequences 
with reddening Ay = 0.25 mag at 2, 8.5, and 25 kpc (as labeled). 
One sees that most stars are at about 10 kpc distance, implying 
that their photometric parallaxes are small, about 0.1 mas. 

empirically from variations in background regions, and 
G = 1 e~ DN _1 is the gain for pipe- line processed HRC 
images). We fit simultaneously for the position and am- 
plitude of each star, but fix the sky to the header value 
mdrizsky, which is determined by the pipeline drizzle 
process (and is very similar to the median or mode of 
all pixel values) . Finally, we correct the po sitions for 
distortion in the HRC using the solution from IAK04L 

To estimate uncertainties on the positions, we used the 
standard technique for \ 2 fitting in which we determined 
the x and y offsets at which x 2 increased by 1. It is 
worth noting that, generally, this is only valid if the fit 
is acceptable, i.e., if x 2 roughly equals the number of 
degrees of freedom. For the brightest stars, we often 
find much larger \ 2 , since the PSF varies due to focus 
changes, etc. One could rescale the errors, but this would 
lead one to overestimate the uncertainties of positions of 
stars relative to each other (which is what enters our 
analysis), since PSF changes affect all stars in the same 
way. Therefore, we retained the formal uncertainties, 
verif ying that they lead to acceptable fits (see below and 
App.[EL2}. 

With positions and uncertainties in hand, we solve 
for the best-fit transformation between the first and the 
other seven exposures for each epoch (for details, see 
App. [X]) , and determine one set of average, distortion- 
corrected positions per epoch. Given our weighting, the 
four brightest stars (A, B, 118, and 1 19) dom inate the 
fits. For the transformation, we follow |AKQl and use a 
6-parameter, bi-linear transformation, as we do for the 
transformation between epochs (see §[4j. This fits for a 
central position, the overall plate-scale and rotation of 
the image, and a second plate-scale and rotation that re- 
flect differences in scale and orientation between the two 
axes (we test variations on this in App. IB.4j) . 

In general, when we included all measurements, the fit 
quality for some transformations was very poor. Inspect- 
ing the outliers, we found that within the fitted region of 
most there were some pixels flagged as bad. These pixels 
are not included in the ePSF fit, so to first order they 
should not affect the result. There is a second-order ef- 



fect, however: that arises because our model ePSF is not 
a perfect match to the real one. As stated above, gener- 
ally, this does not matter for relative positions, since the 
ePSF is wrong in the same way for all stars. But this re- 
lies on the same part of the ePSF being used, which is not 
possible if bad pixels are present. Because of this, we de- 
cided to exclude all position measurements that had bad 
pixels. In addition, we excluded all measurements of star 
114 in epoch 4, where it is close to the HRC occulting 
finger and clearly deviant, and three further significant 
outliers without obvious causes (star 115 in epoch 2, ex- 
posure 2; 118 in 3, 5; and A in 7, 7). 

With the above exclusions, the fits are fairly good, with 
the overall x 2 e( j ranging from 1.00 for epoch 8 to 1.48 for 
epoch 5. Furthermore, all individual sources were found 
to fit well too (for details, including simulated data, see 
App. El and Fig. p. 

2.3. Absolute Astrometry 

While the HRC observations provide very accurate rel- 
ative positions, the precision with which these can be tied 
to the International Celestial Reference System (ICRS) 
is not as good, as it is determined by the less accurately 
known positions of the guide stars used. To improve the 
precision, we use ground-based observa tions to tie our 
meas urements to the UCAC2 catalog (Zacharia s et al.l 
I2004T ) , which is currently the most precise representation 
of the ICRS for stars fainter than those measured by 
Hipparcos. 

We first tied a short FORS image taken on 31 Decem- 
ber 2002, which is very close in time to our third HRC ob- 
servation, to the UCAC2 catalog (corrected to the same 
epoch using the UCAC2 proper motions). We measured 
centroids of 13 UCAC2 stars present on the image and 
corrected for radial distortion using the quadr atic rela- 
tion given bv lJehin. O'Brien, fc SzeiferiJ (|2006f) . A four- 
parameter transformation sufficed, with the offset accu- 
rate to 10 mas, and scale and rotation different from nom- 
inal by -0.028±0.012% and -0?076±0?008, respectively. 
The fit was good, with the residuals consistent with the 
UCAC2 uncertainties. 

We then used 41 stars to transfer the tie to a deeper im- 
age, consisting of twenty 620-s images from 29, 30, and 31 
December 2002. This fit was again good, and the uncer- 
tainties negligible compared to the first step. Finally, we 
tied our astrometry to the HRC using 18 stars for which 
we could obtain accurate FORS positions. Using HRC 
positions evaluated at the FORS epoch using the proper 
motions derived in £|H we find an excellent tie, with off- 
sets accurate to ~ 2 mas and a scale and position angle 
different from the nominal values by —0.034 ± 0.015% 
and — 0?018±0?010, respectively (here, the uncertainties 
take into account that the HRC shows deviations from 
equality of scales and orthogonality; see § [4j) . 

The overall precision with which our measurements are 
on the UCAC2 system is limited by the first step, and is 
~11 mas, simi lar to the precision w ith which UCAC2 is 
on the ICRS (jZacharias et al.ll2004T ). For the scale and 
position angle, the first and last step contribute; we esti- 
mate that the scale is accurate to 0.02% and the position 
angle to 0?011. 

Below, we will give positions relative to star A, fix- 
ing its position to a J2 ooo = 07 h 20 m 24'?4837, <5 J2 ooo = 
— 31°25'51'/786. We will also correct the nominal scale 
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"phot < mas > 

Fig. 4. — Photometric versus astrometric parallaxes for the stars 
from Table [2] for which the astrometric measurement error was less 
than 0.3 mas; the dashed line indicates 7r as t = ""phot- In measuring 
the astrometric parallaxes, we held all reference stars but the one 
labeled fixed at their photometric parallaxes, and then solved for 
the parallax of the reference star in question. 

and position angles of epochs 3 and 7 using the values 
found above. 

2.4. Photometry 

For the HRC photometry, we measured fluxes in the 
averaged flat-fielded images for each epoch (see § I2.2p . 
Since the HRC flux for RX J0720.4-3125 will help con- 
strain the different components of the optical emission 
(|Kaplan et all l2003h . we took care to correct for all 
known systematic e ffects, following the prescription of 
ISirianni et all ((20051) for flat-fielded, non-drizzled images. 
Briefly, we first measured co unts in a range of apertures 
using DA0PH0T (|Stetsonlll987l ). and sky in a 50 to 60 pixel 
annulus (r733-l'/59). We chose an aperture of 3 pixels 
radius (0'.'08) as a compromise between adequate signal- 
to-noise ratio for fainter targets and minimal effects from 
PSF variations over the detector. Next, we corrected 
for pixel-are a variations and charge transfer inefficien- 
cies following iPavlovskv et a l. (2006). We used brighter, 
isolated stars to determine the aperture correction to an 
18-pixel (0'.'48) aperture, and a dded a fixed —0.103 mag 
to correct to "nominal infinity" (|Sirianni et al.ll2005l ) ; the 
latter includes the small, —0.003 mag effect from contam- 
ination of the sky by starlight. Finally, we added the zero 
point of 25.623 to obtain magnitudes in the Vega system. 

We checked for, but found no significant variations 
for any star (for star 113, we rejected epochs 1 and 3, 
since a trail from star A passed over the image). For 
RX J0720.4— 3125, where we list the individual measure- 
ments in Table [TJ the standard deviation is 0.10 mag, 
consistent with the error of 0.09 mag and with numbers 
for other faint sources (e.g., 0.09 and 0.10 mag for stars 
112 and 108, respectively). We list the average magni- 
tudes in Table H 

For STIS, we again used DADPHDT to measure 
fluxes, choosing a 2-pixel (O'.'l) aperture. We 
corrected for charge transfer ineffic iencies following, 
Goudfrooii fe Kimble! {2002, see also iGoudfrooii et al.l 
2006|) . and used brighter stars to correct to a 10 pixel 



(0'.'5) radius aperture. To place our magnitudes roughly 
on the STScI system, we added an additional —0.100 mag 
to "nominal infinity," and used the zero point flux of 



0.9987 x 10~ 19 erg cm" 2 s" 1 



A for a count rate of 



Is" 1 , and the magnitude offset of 21.1, as described in 
the HST data handbook. We stress that while the re- 
sults are fine for our purpose of inferring stellar colors 
(since calibration errors will be taken out), care has to 
be taken in interpreting them in terms of absolute fluxes. 
In particular, for the flux of RX J0720.4-3125 itself, see 

the detailed discussion in lKaplan et alj (120031). 

For the LRIS photometry. iKulkarni fc van Kerkwiikl 
(1998) had found that scattered light from much brighter 
stars made photometry difficult, and they adopted a sim- 
plified PSF fitting technique. In order to measure not 
just isolated stars, but also ones near others, we modified 
their procedure slightly, using the positional information 
from the HRC and picking a moffat function instead of a 
Gaussian as a PSF model. Th e results, listed in Tabled 
are co nsistent with those of IKulkarni fc van Kerkwijkl 
(1998). They are also ro ughly consistent (bright er by 
~ 0.03 mag) with those of IMotch fc Haberli (|1998l ) and 
iMotch et alJ ([20031 ). 

3. PHOTOMETRIC PARALLAXES OF 
BACKGROUND SOURCES 

In measuring the parallax of RX J0720.4— 3125, we 
need to worry about bias by parallaxes of our background 
sources. Typical distances may range from 1 to lOkpc, 
inducing an effect at the 1 to 0.1 mas level. Generally, one 
might expect more nearby sources to bias the fit most, 
since these would be brighter and hence have heavier 
weight. Given this, it is good to have distance estimates, 
which, even if wrong for individual sources, reduce any 
systematic bias. For this purpose, we use our photome- 
try: we estimate temperature from B — R and flux from 
R, and then infer a distance assuming the stars are on 
the main sequence. 

For many stars, we have B and R magnitudes from our 
Keck images f ^2.4|) . Some stars, however, were too faint 
(in particular in B) or too close to brighter stars for reli- 
able photometry. For those, we infer B — R and R from 
the HRC and STIS photometry, using empirical color- 
color relations determined from stars with uncertainties 
in B— R smaller than 0.1 mag; see Fig.[2]and Tabled We 
find that all relations are tight, with root-mean-square 
residuals of less than 0.05 mag. Between 777,475 — R and 
B — R, the relation is almost linear, while the relations 
between 771475 — m 50 and B — R, and 771475 — m 50 and 
m 50 — R are well-described by quadratic functions. By 
way of verification, we compared our result for 771475 — R 
and B — R with the relation exp ected from syntheti c pho- 
tometry on model atmospheres (Siri anni et al.|[2005l ). As 
can be seen in Fig. [2l the agreement is good. 

We proceeded by estimating B—R values from 771475— R 
and 777475 — 77150 for all stars; for a given star, we adopt 
the B — R estimate with the smallest uncertainty (where 
the uncertainty includes in quadrature the measurement 
uncertainty and the scatter around the required transfor- 
mation). For objects for which B — R was inferred from 
771475 — to 50j we estimated R from 77750 using the relation 
between 777475 — 77750 and 77750 — R. 

In Fig. [3j we show our adopted B — R and R estimates 
in a color-magnitude diagram. Also shown are expected 
colors and magnitudes for main-sequence stars at various 
distances, for a reddening of Ay — 0.25 (see below). For 
thes e, we used the absolute magnitudes My and colors 
from lCoxl (|2000l : specifically, we obtained My and B — V 
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from Table 15.7, and converted to Mr and B — R us- 
ing Table 15.11; we cannot use Table 15.7 directly since 
our R-band magnitudes are in the Kron-Cousins system) . 
One sees that if our objects are main-sequence stars, they 
lie between 2 and 20 kpc; the estimated photometric par- 
allaxes are listed in Table [2] We note that it is very un- 
likely that they are not main-sequence stars: if they were 
giants, they would be well outside the Galaxy, while if 
they were white dwarfs, we would find a significant par- 
allax. 

The photometric parallaxes for the brighter stars, 
which carry most weight in our astrometry, are all around 
0.1 mas, and thus even large fractional errors in these 
will not influence our conclusion much. Nevertheless, 
it is worth briefly considering the main sources of er- 
ror. Likely these are the assumed reddening and metal- 
licity. F or the reddening we us ed Eb-v = 0.08, de- 
rived by iMotch fe Haberll (| 19981) from UBV photom- 
etry of stars near RX J0720.4-3125. This is con- 
sistent with the Eb~v — 0-07 inferred from the ob- 
served B — R colors of stars A and B and their G2 
and G5 spectral ty pes (as i nferred from the spect ra of 
lHaberl et all Il997t see also IMotch fe Haberll [l998h . It 
is somewhat lower than the total line-of-sight extinction 
Eb- v = 0.14 estimated from background infrared emis- 
sion (Schlcg el. Finkbeiner. fc D avis 199§|). 

The above estimates assume solar metallicity. A sub- 
solar metallicity may be more likely, however, since 
the background sources are relatively high above the 
plane and at large Galactocentric radius: for (I, b) — 
(244?16, -8?16) and at 10 kpc, one finds z = 1.4 
above the Galactic plane and Galactocentric radius w — 
16 kpc, (we ignored tha t the plane is warpe d down 
by - 3° at t = 240°; iMomanv et~aH I2006D . This 
may also influence the reddening: for [Fe/H] = —1.0, 
IMotch fc Haberll (|1998f ) infer a lower reddening Eb-v — 
0.04 from their UBV photometry, while from stars A and 
B one infers a higher reddening Eb-v — 0.14 (since they 
would be intrinsically bluer). 

To estimate the maximum effect of sub-solar metal- 
licity, we combined the increased intrinsic blueness, de- 
creased luminosity, and increased reddening for star A. 
We found that a change of 0.6 mag in distance modulus, 
or a factor 1.3 in photometric parallax. For stars of later 
spectral type, the effect is less. Since the brighter back- 
ground stars are so distant, the additional uncertainty 
in the parallax of RX J0720.4— 3125 should be less than 
0.1 mas. 

4. MEASURING THE PARALLAX 

Our goal is to measure the proper motion and parallax 
of RX J0720.4-3125, but this has to be done relative 
to other stars, for which the positions, proper motions 
and parallaxes are not known a priori, and on images 
for which our knowledge of the positions, rotations, and 
scales is insufficiently accurate. In our procedure, we fit 
simultaneously for our target parameters and all other 
parameters. We have to fix some of those, however, since 
otherwise the fit is degenerate (e.g., a net shift between 
epochs cannot be distinguished from a proper motion 
component common to all stars). We made the following 
choices. 

First, we chose epoch 3 to have a known plate-scale and 
position angle, with the offsets from the header values 
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Fig. 5. — Positions of the reference stars at every measured epoch, 
with the predicted track based on the best-fit proper motion and 
photometric parallax indicated by the solid line. The diamond 
indicates the fitted position at the epoch of the first observation. 

(Table [1} set by our ground-based absolute astrometry 
(§ 12. 3p . This sets the absolute scale and orientation of 
the data. 

Second, we also chose epoch 7 to have known plate- 
scale and position angle, using the header values with 
the same offsets as for epoch 3. This ensures that the 
scale and position angle do not drift with time (i.e., it 
avoids net expansion or net rotation) . We chose this pair 
of epochs since they are at almost the same parallactic 
angle, so their tie cannot influence the parallax. The 
tie will influence the orientation and scale of the proper 
motions, leading to an uncertainty of 0?05 and 0.01% (as 
estimated from the differences in orientation and scale 
found for other epochs below). This effect, however, is 
smaller than that of our next step. 

Third, we fix the position of star A to that determined 
from the absolute astrometry, and set its proper motion 
to 0. This sets the reference position and fixes the net 
proper motion. Given the proper motions of other dis- 
tant stars, we estimate that it introduces an uncertainty 
in our net proper motion of ~ 1 masyr -1 in each coordi- 
nate. This is very small compared to the proper motion 
of RX J0720.4— 3125, but comparable to the formal un- 
certainty, and thus should be taken into account. We 
note that the choice of star A is arbitrary, and one could 
have picked another star or an ensemble. An ensem- 
ble might have similar proper motion, however, and our 
choice has the advantage of being simple, and of star A 
being a good reference object in that it is bright, has a 
small photometric parallax, and is located near the cen- 
ter of the images. 

Finally, to fix the mean parallax (which would be indis- 
tinguishable from an epoch-dependent shift in the point- 
ing), we fix the parallaxes of all stars but our target to the 
photometric parallaxes from Table [3J To verify whether 
these were reliable, we computed solutions where we fit- 
ted for the parallax not only of our target but also of one 
more star. Cycling this additional star among all refer- 
ence stars, we find good agreement between the photo- 
metric and astrometric parallaxes (see Fig.|4|. Even the 
most deviant points is only 2.3<r away, and equality has 
X 2 = 22.5 for 22 degrees of freedom (22 measurements, 
no free parameters). Thus, our photometric parallaxes 
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TABLE 4 

ASTROMETRY OF REFERENCE SOURCES 







^instr 


Act 


A<5 


Pa 






ID 


N u 

epoch. 


(mag) 


(arcsec) 


(arcsec) 


(mas yr -1 ) 


(mas yr~ 1 ) 




X 


8 


-5.7 


6.1219(8) 


1.7026(9) 


-93.9(12) 


52.8(13) 


1.13 


A 


8 


-12.0 














0.67 


13 


8 


-11.8 


9.46094(12) 


-8.01896(11) 


-0.45(17) 


-0.18(17) 


0.50 


104 


7 


-9.2 


13.92144(16) 


-2.65145(15) 


0.4(2) 


0.8(2) 


0.44 


105 


8 


-6.4 


4.2079(5) 


4.0488(5) 


-0.2(7) 


-0.7(7) 


1.83 


106 


8 


-6.5 


3.4296(5) 


5.5332(5) 


2.9(7) 


-1.8(7) 


0.93 


107 


8 


-9.3 


2.03280(16) 


10.34588(16) 


1.7(3) 


-0.3(3) 


0.51 


108 


7 


-5.4 


3.8609(14) 


11.8761(13) 


-0.2(19) 


-0.9(18) 


1.19 


109 


3 


-7.8 


6.1437(5) 


15.4211(5) 


1.3(7) 


-1.2(7) 


1.39 


110 


5 


-10.2 


-2.9762(3) 


13.0787(3) 


-0.2(4) 


6.3(4) 


0.53 


111 


7 


—9.0 






—0.3(3) 


0.3(3) 


0.98 


112 


8 


-5.8 


-1.6390(8) 


1.4360(8) 


-0.3(12) 


0.2(13) 


0.96 


113 


8 


-5.5 


0.1059(10) 


-1.3660(9) 


2.3(16) 


0.9(15) 


0.46 


114 


7 


-10.2 


-5.84398(9) 


0.22906(9) 


-0.25(14) 


1.28(14) 


0.99 


115 


8 


-10.9 


-6.88576(9) 


-2.06279(8) 


-0.69(13) 


0.83(13) 


0.62 


116 


6 


-9.8 


-10.62105(13) 


-1.09834(13) 


-1.8(2) 


-0.0(2) 


0.12 


118 


8 


-12.0 


-6.52202(9) 


-8.86165(9) 


-0.68(15) 


-0.86(15) 


0.60 


119 


1 


-12.8 


-10.16052(10) 


-10.01889(10) 


-0.86(19) 


0.16(19) 


0.27 


120 


1 


-8.3 


-2.8154(4) 


-12.5970(4) 


-4.7(5) 


-6.1(5) 


0.95 


121 


1 


-7.2 


-1.6247(5) 


-13.9458(6) 


-0.4(7) 


0.4(7) 


1.81 


122 


1 


-8.8 


-6.5404(4) 


-15.8855(4) 


-2.5(5) 


1.0(5) 


0.19 


123 


8 


-7.0 


11.9215(4) 


-7.3533(4) 


-0.6(5) 


0.3(5) 


2.01 


121 


8 


-5.9 


12.6335(8) 


-7.6056(8) 


-0.5(11) 


1.3(11) 


1.00 



Note. — Columns are: ID: star number, with letters following Habcrl et al. (1997). Source 
X is RX J0720.4— 3125; for its parallax, see Table [6] N cpoc ^: number of epochs for which a 



source could be measured. m,j nstr : instrumental magnitudes defined as — 2.51og 10 a, where a 
is the average amplitude of the ePSF required to fit a giv en star in a single exposure. Act, A5: 
position offsets at MJD 52645.1, relative to star A (§ 12,31 1. excluding parallactic offsets. fi a , fig: 
proper motion relative to that of star A. The systematic uncertainty on the proper motion, 
from assuming that star A has fi = 0, is ~ 1 mas yr — 1 . X 2 c d : r °duced X 2 values for fitting, 
with degrees of freedom 2N epoc h — jV par , where the number of parameters 7V par = 5 for 
RX J0720.4— 3125, for A, and 4 for all other stars (note that these X 2 e d values for individual 
sources are not robust estimates; see App. [X). All uncertainties are at the ler level; entries 
with no uncertainties were fixed at the given values. 



TABLE 5 

Parameters of the Epoch Transformations 







x 


yo 




A 


4>y 


Ay 




Epoch 




(pix) 


(pix) 


(10~ 3 deg) 


(io- 5 ) 


(10~ 3 deg) 


(io- 5 ) 


•^red 


1 


17 


556.164(3) 


496.457(3) 


99.3(9) 


-30.9(12) 


-2.9(11) 


-1(2) 


0.30 


2 


21 


521.164(3) 


444.095(3) 


90.5(10) 


-33.6(19) 


8.0(15) 


-9(3) 


0.35 


3 


18 


461.598(3) 


517.035(2) 


17.71 


-34 








0.95 


1 


21 


491.148(3) 


573.595(3) 


12.7(9) 


-21.0(16) 


11.7(13) 


-17(2) 


0.55 


5 


17 


581.178(3) 


492.582(3) 


9.1(10) 


-29.4(11) 


1.0(12) 


-6(2) 


1.63 


6 


20 


546.786(3) 


471.803(3) 


14.6(6) 


-28.1(12) 


14.7(9) 


-20.1(16) 


0.85 


7 


19 


467.461(3) 


520.372(3) 


17.71 


-34 








0.66 


8 


21 


509.103(4) 


581.869(4) 


10.0(8) 


-20.1(16) 


16.4(12) 


-19(2) 


0.58 



Note. — Columns are: Epoch: epoch number. 7V s t ar : number of stars included. xo,j/o : reference 
position, which is the model position of star A in dedistorted coordinates. <fi: deviation of the best-fit 
position angle relative to the nominal orientation listed in Table [T] A: deviation of the best-fit scale 
relative to the nominal plate scale of 28.27 mas pixel -1 . (p y ,X y : deviations from orthogonality and 
equality of scale; see Eq. lA5l X 2 c d : reduced x 2 values for fitting, with degrees of freedom 2iVst a r — N paI , 
where the number of parameters JV par = 2 for epochs 3 and 7, and 6 for all others (note that these 
X 2 ed values for individual epochs are not robust estimates; see App. |XJ. All uncertainties are at the 
la level; entries with n o un certainties were fixed at the given values, which were determined from our 
absolute astrometry (§ 12.30 . 
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TABLE 6 
Motion of RX J0720.4-3125 



Parameter Best-fit Values p w ^ 



a J20 oo b 07 h 20 m 24 ! :9620 ± 0?0009 -0.08 

<5j2ooo b -31°25'50'.'083±0'.'011 -0.10 

/i„ (mas jt- 1 ) -93.9 ±1.2 0.20 

US (mas yr" 1 ) 52.8 ± 1.3 0.25 

7T (mas) 2.77 ± 0.89 ■■■ 

D (pc) 361^ 

fi (mas yr" 1 ) . 107.8 ± 1.2 

PA (deg) -60.6 ± 0.8 

v ± (kms- 1 ).. 185±H 



Note. — Errors for a an d 8 are determined by the 
absolute astrometry (§ 12. 31 1. Errors for fj, a , fig, and 
7r are formal uncertainties at the lcr/68% confidence 
level. There may be additional syste matic uncertain- 
ties, of at most 0.4 mas for it (see App. [B~6t and about 
1 mas yr — 1 for fj, a and fig (see §|4j. Best-fit values and 
errors for the other parameters are derived from those 
for /i a , fig and it, ignoring the possible systematic un- 
certainties. 

a The correlation coefficient between each parameter 
and the parallax n. 

b At epoch MJD 52645.1 and for equinox J2000, ex- 
cluding any parallactic offset. 



appear reliable, at least on average. 

With all parallaxes fixed except for that of 
RX J0720.4-3125, we derived the fit given in TablesHEl 
and [SI The quality of the fit was good, with x 2 — 192.1 
for 179 degrees of freedom (306 data points and 129 pa- 
rameters), and most stars fit quite well, as one can see 
in Fig. 03 where we compare the measured positions for 
the reference stars with our fit. Examining each epoch 
independently (Table [5]), we find that only epoch 5 has 
Xred significantly greater than unity. Interestingly, this 
epoch also had the poorest \ 2 m the registration (§ I2.2p . 
Looking in detail, we found that the ePSF fits generally 
had worse x 2 , and that the small-aperture instrumental 
magnitudes were fainter by ^0.2 mag compared to those 
of the other epochs (this did not affect our photometry, 
since the aperture corrections compensated for it). Vi- 
sually, the PSF in the images from epoch 5 appears less 
sharp, and thus we believe that this epoch suffered from 
larger than average focus variations. We note, though, 
that while this explains the fainter small-aperture mag- 
nitudes and poor ePSF fits, it is unclear why the registra- 
tion or astrometry should be significantly poorer. Over- 
all, we felt there was insufficient reason to reject epoch 5 
and it is therefore included in our final analysis (we ex- 
amine what happens if we exclude epoch 5 — and other 
epochs — in App. IB. 5]) . 

For RX J0720.4-3125, we find a parallax tt = 2.77 ± 
0.89 mas; the full set of fitted parameters is listed in Ta- 
ble[Hl and the fit is shown in Fig. [SI The fit is good, with 
Ared = 1-13 (for 11 degrees of freedom; 16 measurements 
and 5 parameters). From Fig. one sees that also for 
RX J0720.4— 3125 itself epoch 5 fits worse than any of 
the others, with a net deviation of 2.4cr, while the next 
worst is epoch 8 (1.6<r). We return to this in App. [Bl 
where we try to estimate systematic errors in the paral- 
lax; we find that these are at most 0.4 mas, i.e., less than 
half the measurement error. 
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Fig. 6.— Parallax of RX J0720.4-3125. Upper left: parallactic 
ellipse showing motion in Right Ascension versus that in Declina- 
tion with the proper motion removed. The filled circles are the 
predicted positions of RX J0720.4— 3125, and the open diamonds 
are the data. The data are connected to the predicted positions at 
the same epoch. Upper right: motion in Declination versus time 
with the proper motion removed. Lower left: motion in Right As- 
cension versus time with the proper motion removed. Lower right: 
motion in Right Ascension and Declination including the proper 
motion. 



5. DISCUSSION 

The distance d = 360tgQ° pc we measure is roughly 
consistent wi t h the estimate of 250 ± 25 pc made by 
iPosselt et all (2006) comparing the hydrogen column 
density to RX J0720.4-3125 — as inferred from fits 
to its X-ray spectrum — with the run of hydrogen col- 
umn with distance — as inferred from sodium columns 
to stars with known distances. Conversely, our parallax 
measurement indicates that the column density inferred 
from the spectral fits, which necessarily depends on the 
form assumed for the intrinsic spectrum, is reasonable 
(if perhaps slightly low; a distance slightly closer than 
the parallax distanc e is also found for RX J1856.5— 3754: 
IPosselt et all (|2006l ) find d = 135 ± 25 pc, while our pre- 
liminary parallax yields 167^15 P c )- 

Our parallax implies that RX J0720.4— 3125 is roughly 
a factor two more distant than RX J18 56.5— 3754 . This 
agrees well with the simple estimate of lKvKA02t which 
was based on the first-order assumptions that the optical 
flux for different sources scaled as F v oc T(R/d) 2 , that 
the radii were similar, and that the temperature T in 
the region emitting optical photons scaled with the tem- 
perature determined from fits to the X-ray spectra. Al- 
though this agreement may be a coincidence, it suggests 
that the mismatch between black-body fits to the X-ray 
and optical emission arises because the surface does not 
emit like a black body (i.e., that the temperatures in 
the layers of the atmosphere emitting X-ray and optical 
photons are different but rela ted; e.g.. lMotch et al.|[2003l : 
iZane. Turolla. fc Drake|[2004T) . and not because the opti- 
cal and X-ra y emission originate in s eparate hot and cool 
regions (e.g. JBraie fc Romanil [20021 ). 

While the X-ray spectrum and flux of 
RX J0720. 4-3125 had been observed to be rela- 
tively constant for years after its discovery (like for 
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Fig. 7. — Origins of two nearby isolated neutron stars. This figure, adapted from Fig. 29 of de Zee uw et al.l JT999), shows a view looking 
down on the Galactic plane (right panel, with Galactic longitude I as indicated) and a view acro ss the plane (left pan el, with Galactic 
latitude b as indicated). The filled circles show kinematically detect ed OB assoc iations (filled circles; dc Zeeuw et al. 1999) alo ng with their 
stream ing motions, and the scattered dots represent the Gould Belt (Olano 1982). The present locations of RX J1856.5— 3754 IIKaplan et al.l 
2002b; van Kerkwijk et al. 2006, in prep.) and RX J0720.4— 3125 (this work) are indicated, as are trajectories back in time for the no mina l 
distance and proper motion, assuming radial velocity ti rat j = 0, as well as using the lcr distance limits and v rat j = ±0.935tix (see ij 15.11 1. 
The tick marks along the trajectories occur every 0.5 Myr; Galactic acceleration is negligible over these timescales. The scale for velocity 
of 20 km s _1 for the associations is shown along the bottom, while the linear scale of 500 pc is along the right edge. For reference, a source 
with a velocity of 20 km s — 1 will traverse 500 pc in 25 Myr. 



the other INS), observations in 2 002 and 2003 s howed 
a surprising sp ectral change (|de Vries et alj l2004t 
IVink et aTll2004h . The spectrum hardened significantly, 
although the flux stayed relatively constant. Our ACS 
observations occurred during the same time span. In 
§ 12.41 we looked for variability in the optical flux of 
RX J0720.4-3125, but did not find any. This may 
not be entirely unexpected, however, if the optical flux 
indeed results from the Rayleigh-Jeans tail of some 
hot region (and scales as R 2 T, as above). Over the 
course of our observations, the temperature of the X-ray 
blackbody changed by 6%, from kT = 88.3 ± 0.3 keV 
(2002 November) to 93.8 ± 0.4 keV (2004 May), while 
the blackbody an gular diameter dro pped ~ 6% over 
the same period (Ha berl et alj 1200(1 ). This implies a 
decrease in optical flux of ~ 6% = 0.07 mag, which is 
smaller than our photometric uncertainty of 0.09 mag 
(Tabic QJ. Even binning the observations we do not see 
significant variability: comparing the mean of epochs 
1-4 with 5-8 there is a decrease of 0.08 mag, but the 
uncertainties on each mean are 0.05 mag. 

5.1. Origin and Age 



With the distance and proper motion of 
RX J0720.4— 3125, we can estimate the space mo- 
tion and try to determine its origin and age. As the 
simplest estimate, we determine the time required for 
RX J0720.4— 3125 to reach its current location out of 
the Galactic plane. The current Galactic latitude is 
b = —8.16°, and the proper motion in that coordinate is 
/!/, = —59.3 mas yr _1 . So, assuming it started at b = 0, 
it has taken ~ b/fii, — 5 x 10 5 yr to reach that height. 
However, this value is rather inaccurate. At a distance 
of 360 pc, RX J0720.4-3125 is at z - z Q = -50.7 pc, 
where z is the height above the Galactic plane, an d the 
Sun is at z^, sa 10 pc c ompared to local OB stars (jReedl 
ri997tlElias et~all[200l . So RX J0720.4-3125 is actually 
at z w —40 pc below the mid-plane. Given that this is 
comp arable to th e local scale-heigh t of 30-40 pc of OB 
stars (Reed 2000; Elias et al.ll2006f ). we cannot actually 
use this method for a useful age estimate, but can only 
derive a rough upper limit of ~ 1 Myr. 

We can derive a more meaningful estimate from con- 
sidering not simply the scale height of the OB pop- 
ulation, but instead the actual locations of individual 
associations. In Fig. [71 we show the nearby Galac- 
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tic OB associations (taken from Ide Zeeuw et al.l [19991. 
and overdraw the current location and the previous 
motion of RX J0720.4— 3125 for the distance corre- 
sponding to the best-fit parallax and zero radial veloc- 
ity «rad- We also show the result for distances corre- 
sponding to 7r ± ov, and for t> rac i = ±0.935u^, where 
v± = 185 kms" 1 is the velocity in the plane of the sky 
inferred from the proper motion and best-fit parallax, 
and the numerical factor cot(cos _1 0.683) = 0.934 cor- 
responds to the expected la range in v rac i for random 
orientations. In our estimates of the space velocity, we 
corrected for the Sun's motion relative to the Local Stan- 
dard of Rest, of (U V, W ) e = (10.00,5.25,7.17) kms" 1 
(|Dehnen fc Binnevl 119981 ). as well as for differential 
galactic rotation, using P ort constants (A, B) = 
(14.82,-12.37) kms-ikpc" 1 (jFeast fc Whitelocld[l99"l . 
For comparison, we also show the location and 
motion of RX J1856.5— 3754 (using [n ai fis] = 
[0.327, -0.059] mas yr" 1 and tt = 6.0 ± 0.6 mas [van 
Kerkwijk, Kaplan, & Anderson 2007, in preparation!). 

From Fig. one sees that, as discussed in IWalterl 
(|2001h . RX J1856. 5-3754 plausibly came from the Up- 
per Scorpius or Upper Scorpius Lupus OB association 
~ 0.4 Myr ago For RX J0720.4-3125 as suggested 
bv iMotch etahl (|2003l ) and [Kaplan! ([20041 ) . an origin in 
the Trumpler 10 (Tr 10) OB association ~ 0.7 Myr ago 
seems likely. This is plausible in £, b, and distance, with 
little radial velocity required (values between —20 and 
+50 km s" 1 are compatible with the nominal distance). 
An origin in the Vela OB2 association is less likely but 
not impossible; it would imply a slightly smaller age 
(~ 0.6 Myr). For a substantially larger radial velocity 
(> 500 km s~ ), origins in the Lower Centaurus Crux, 
Upper Centaurus Lupus, or Upper Scorpius OB asso- 
ciation are possible. Even for those, however, the age 
would still be below 1 Myr (as it has to be from the 
the proper motion in galactic latitude, see above; Fig. [7] 
shows that an age in excess of 1 Myr is excluded also for 
RX J1856. 5-3754). 

Overall, among all of the OB associations, an origin 
in Tr 10 seems most likely, with RX J0720.4-3125 ap- 
proaching within ~ 17 pc of the core (the radius of the 
association is ~ 30 pc). This should be compared to 
70 pc for Vela OB2, although this is a larger association 
with a radius of ~ 55 pc. Of course, the fact that the 
trajectory points back to Tr 10 could be a coincidence. 
To estimate the probability of such a coincidence, we 
did a simulation where we assumed proper motions with 
the observed magnitude but with random direction on 
the sky. For each direction, we computed the minimum 
possi ble separation a c hieved among all OB associations 
from Ide Zeeuw et"aT1 (|1999() . as well as the required ra- 
dial velocity. We find that we can get RX J0720.4-3125 
to within 25 pc (comparable to the closest approach to 
Tr 10) of any OB association in only 5% of our trials, 
and in only 2% of our trials if we restrict the radial 
velocity to |i> ra d| < 350 km s _1 (which is almost 2v±). 
If we allow for the fact that some OB associations are 
larger than others, we find that in only 10% of our tri- 
als does RX J0720.4— 3125 approach within the nominal 
radius, and in only 5% of our trials for the restricted 
radial-velocity range. If we vary also the magnitude of 
the proper motion, we obtain similar results. Therefore, 



given the small size of Tr 10 and the closeness of the ap- 
proach, we consider it likely that RX J0720.4— 3125 came 
from Tr 10. 

If we accept an origin in Tr 10, the main uncertainty 
in the age comes from our distance uncertainty. To 
assess this, we drew parallaxes for RX J0720.4— 3125 
from a normal distribution corresponding to our mea- 
surement uncertainty, and determined what radial ve- 
locity and age gave the closest approach to Tr 10. We 
found, as expected, a most probable age of 0.7-0.8 Myr, 
and a linear relation between parallax and age 4 such that 
Age = 0.7 (ir/2.8 mas) Myr. We thus infer an la range 
in age of 0.5-1.0 Myr. Larger parallaxes imply more neg- 
ative radial velocities (approaching —100 km s _1 ), while 
small parallaxes require high positive radial velocities (up 
to 400 km s _1 for tt < 2 mas) to get RX J0720.4-3125 
to the distance of Tr 10 in the short time allowed. 

The age we infer for RX J0720. 4-3125 is compa- 
rable to simple estimat es based on models of neu- 
tron star cooling (e.g., IKaplan et al.1 l2002al ). but it 
poses a pair of puzzles. The first is that the spin- 
down age for RX J0720.4-3125 of 1.9 Myr is a fac- 
tor of 2-3 larger t han either the cooling age or our 
kinematic age (see IKaplan fc van Kerkwiikl l2005al and 
van Kerkwijk fc Kaplanl^OOlT who also discuss possible 
explanations). In this respect, RX J0720.4— 3125 may 
not be unique: the hotter neutron star RX J1308. 6+2127 
has a similarly long spin-do wn age of 1.5 Myr 
(Kapl an fe van Kerkwijk 2005b). The second possible 
puzzle is that, kinematically, RX J0720.4— 3125 ap- 
pears older than RX J1856.5— 3754, while its tempera- 
ture is higher. Thus, the sources likely differ in some 
other aspect. The inferred magnetic field str engths are 
indeed different (~ 2 vs. < 1 x 10 13 G; lKvKA02t 
IKaplan fe van Kerk wijk 2005a), but this may not suf- 
fice. Instead, it may reflect a small difference in mass, 
to which c ooling appears to be especi ally sensitive (for a 
review, see lYakovlev fe PethicldEiol . 

6. CONCLUSIONS 

We have measured the geometric parallax to the neu- 
tron star RX J0720.4-3125 with HST to a precision of 
30%, with an additional systematic uncertainty of ~ 15%. 
To our knowledge, at B = 26.6 mag, RX J0720.4-3125 
is the faintest optical object for which a parallax has 
been measured. While the measurement is too uncertain 
to lead to useful constraints on the radius (which would 
also require better model atmospheres), it helps greatly 
to set the demographic scale of the INS as a population 
relative to that of the radio pulsars, for which parallaxes 
of this magnitude have been measured routinely 5 (using 
Very Long Baseline Interferometry [e.g., IChatteriee et al.l 



20041 for regular pulsar s and "timing" parallaxes [e.g., 



_I E_ 

van Straten et al.l 12001] for millisecond pulsars). Fur- 



thermore, the astrometry shows that the space veloci- 
ties of RX J0720.4-3 125 and RX J1856.5-3754 are typ- 
ical for radio pulsars ([Ar zouma nian. Chernoff. fc Corded 
l27)0ilFaucher^I guere fc Kaspill200o^ 



4 A linear relation is expected since the pulsar needs to traverse 
a certain physical distance between the line of sight and the cluster, 
and, for given proper motion, its space velocity perpendicular to 
the line of sight depends linearly on the distance^ 

5 See ["ittp://www. astro. Cornell. edu~shami/psrvlb/parallax. html http://v 
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The prospect for small improvements of the parallax 
are good, as a few measurements with a large time base- 
line will significantly constrain the proper motions of 
RX J0720.4— 3125 and the reference stars, which will 
improve the epoch registration, reduce the systematic 
uncertainties (App. |B~6|) . and reduce the parallax uncer- 
tainty through its covariance with the proper motion and 
reference position (Table [6]). With such measurement, to- 
tal (statistical plus systematic) uncertainties around 20% 
may be possible. To do much better, however, would re- 
quire further intense observations, and more than 100 
orbits of HST. 

Until a more accurate distance is available, the best 
constraints will probably arise if the data are com- 
bined with our improved distance for RX J1856.5— 3754 
(van Kerkwijk et al. 2006, in prep.). While 
RX J0720.4— 3125 has the poorer distance, in other ways 
the modeling is more constrained, since we know the spin 
period (but see lTiengo fc Mereghettill2~007h . have an esti- 
mate of its dipole magnetic field strength, and can use the 
variation with viewing geometry a nd the broad absorp - 
tion feature in the X-ray spectrum (|Haberl et al. 2004a). 

At present, the most useful aspect of our measurement 
may be that they provide the necessary input for plac- 
ing RX J0720.4— 3125 on a cooling diagram; previously, 



this was done with an erroneous spin-down age, which 
disagrees with our kinematic age and l eads to the objec t 
appearing hotter than predicted fe.g. JPage et al.i r2004) . 
Of course, for proper comparison, we must be careful 
to construct cooling curves that reflect what we know 
about the surface composition (partially -ionized hydro- 



gen? see Ivan Kerkwijk fc Kaplan! 120061) and magn etic 
field (~ 2 x 10 13 G; iKaplan fc van Kerkwijk! l2005al) of 
RX J0720.4-3125. 
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APPENDIX 

A. DETAILS OF THE PARALLAX FITTING 

Our method of determining the parallax from the position measurements in individual exposures was described 
briefly in § |U Here, we describe it in more detail. 

We measured positions using the effective-PSF (ePSF) technique, as discussed in lAK04j 6 and ? 12.21 and this gives us 
64 measurements of up to 23 stars. In order to determine the parallax of RX J0720.4— 3125, it is necessary to transform 
and combine those measurements. We did this in two stages: (i) combine measurements from different exposures in 
each epoch; and (ii) fit the resulting averages with an astrometric model. 

The first stage starts with distortion-corrected positions (x, y)Ees and the associated uncertainties {<T x ,<j y )Ees, where 
E = 1 . . . A^poch is the epoch, e — 1 . . . A oxp is the exposure, and s = 1 . . . N sta , T is the star number, and the uncertainties 
are inferred from the ePSF fit (§ 12. 2ft . For each epoch E, we determine average transformations to a common frame by 
solving for six parameters {(Ace, Ay)E e , ipEe, VEei ipy,Ee, f]y,Ee\-, which relate averaged positions (x, y) e s in the common 
frame to the model positions (x, y)Ees in each exposure, by 

= (1 + TjEe) [{xEs - Ax Be ) COS IpEe + (yEs ~ Aj/ Be ) sin 1p E e] , 
: (1 + TjEe) [—(XEs - As £e ) sin %p E e + (VEs ~ &VEe) COS IpEe] , 

: (1 + 11y,Ee) [-XEes Sfrllpy^e + VEes COS1p y , E e] ■ (Al) 

Here, (Ax, Ay)^ e are the offsets of each exposure compared to the common frame, ipEe is the difference in rotation, 
rjEe is the difference in plate-scale, and ipy.Ee and Vy,Ee are additional parameters that represent the non-orthogonality 
of the x and y axes and the difference between the scales of those axes. Note that while we use all six parameters in 
general, we also made trials with the scale difference and non-orthogonality fixed to zero; see App. IB. 41 

In our fits, we take the common frame to be the first exposure, i.e., for this exposure all transformation parameters 
are zero, and one has [x,y) E is = {p, v)es- We solve for the remaining transformation param eters and the aver age 
positions (x,y)Es at the same time, using a numerical \ 2 minimization routine (mrqmin from lPress et al.lll992D to 
minimize, 

/ X E es - X E es \ , / VEes ~ VEe 



XEe 

VEes 

]jEes 



EE 

s=l e=l 



&x,Ees 



<3y,Ee 



(A2) 



To verify whether our fits our reasonable, we use the global x 2 value, and compute \ 2 values for individual exposures 
and stars. For a more direct comparison with our input uncertainties dEes, we also calculate the standard deviation 
for individual sources, 



„2 

b x.Es 



1 



JVexp 

e=l 



{x E es - X Es ) 



(A3) 



6 For the ePSF mode l, distortion solution, and 
[nttp:/ /spacibm.rice.edu~jay/HRC http:/ /spacibm. rice.edu/~jay/HRC. 



astrometric software described by 1AKQ4. see 
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Fig . 8, — Positional uncertainty in pixels versus instrumental magnitude for the HRC data — compare to Fig. 1 from Kv KA02 and Fig. 2 
from I AK041 We show the standard deviations of the individual positions around the average position for each epoch (Eg. IA3t . in both x 
(crosses) and y (open circles). We also show the root- mean-square average of all the x and y standard deviations (filled squares). Note that 
for the brightest stars, the standard deviations are e xpec ted to underestimate the true measurement uncertainties (App. [Ajl. The curve 
gives the relation expe cted from our simulations (Eg. IBlb . The simulations may also underestimate the true uncertainty at the brightest 
magnitudes (Add. [B72ll . 

and the same for y. Since iV exp is sometimes small for a given star, the standard deviation is not always a good 
estimator of the true uncertainty, but on average we have found that it is quite reasonable (see § IB. 21 and Fig. [5]). 

In considering x 2 values and standard deviations for individual exposures and objects, one has to keep in mind that 
the transformation parameters and the average positions are not independent. In our case, the four brightest stars (A, 
B, 118, and 119) dominate the fit. These provide eight measurements per exposure (4 in each coordinate), but there 
are six free transformation parameters, and hence the fit will partly adjust to remove measurement errors. As a result, 
the x 2 for these sources are often much less than 27V oxp — 2, the value one would expect if the transformation were 
independent of the average positions. For fainter stars, however, which have much less effect on the transformation, 
ones does expect (and we find) x 2 — 2-/V cxp — 2. Of course, for the same reason, the standard deviations calculated 
using Eg. IA3I will underestimate the true measurement uncertainties for the brightest stars. 

Since we fit for all parameters at the same time, the above covariances are taken into account automatically. For 
instance, for star 119, our fit yields errors oes on the average position similar to the input uncertainty OEes for a single 
measurement; this is a consequence of the fact that substantial variation in the average position can be compensated for 
by changes in the transformation parameters. In contrast, for fainter stars, we find oes — &Ees/ \/ N cxp , as one would 
expect for averaging -/V oxp exposures using a fixed transformation. Nevertheless, we worried about our transformation 
being dominated by just a few stars, especially since their formal measurement uncertainties are very small, as low 
as 0.003 pix for star 119, at which level systematic effects may well dominate. Thus, in App. IB. 21 we verify that our 
results are robust to variations in the weight assigned to the brightest stars. 

After the combination of the exposures, we have a set of {x,y)ss along with their uncertainties. We also have a 
set of times tg (which determi ne parallac tic offsets [Aa^, AS^e computed using either the JPL DE200 ephemeris 
or the approximate formulae of Cox 2000, p. 670; both methods gave identical final results) and initial guesses for 
the position angle PA^ (based on the data headers; Table [1} and plate-scale scales (28.27 mas pixel" 1 ; EIqI. In 
our second stage of the parallax determination, we want to solve for the full set of stellar parameters S s and epoch 
transformations parameters £e- 

The stellar parameters S s are the reference celestial position (Aa, A8) s of star s relative to star A at time to (note 
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that we chose equal scales for Aa and AS, so no cos 5 term is needed), the proper motion {^ a ,^s)s, and the parallax 
7r s (for individual stars, the proper motion or parallax may be fixed or set to zero). In terms of these parameters, the 
celestial position of star s at epoch E is given by, 



Aa Es = Aa s + ^ a . s {t E - t ) + Tr s Aa W)E , 
AS Es = AS S + fj,s, s (t E - t ) + ^A^s, 



(A4) 



For the exposures, the parameters £ E are the position (xo,?/o)-E on the detector of the reference position (i.e., 
[Aa, AS] = [0,0], or the [model] position of star A), the difference </>s between the initial guess PAs for the position 
angle and the fitted value, the difference X E between the initial guess scales for the plate-scale and the fitted value, 
and additional parameters 4>y,E and X y<E that represent the non-orthogonality of the x and y axes and the difference 
between the scales of those axes. With these parameters, the transformation from the celestial position of each star 
at each epoch (Ad, AS) Es to the model positions {x,y) Es is given by, 

x' Es =scale^ 1 (l + \ E ) —Aa Es cos (PA B + 4> E ) + AS Es sin (PA E + 4>i 



VEs 



scale/ (1 + As) Aa Es sin (PAs + <p E ) + AS Es cos (PAs + 0s) 



XEs =Xo, E + X Es , 
2/Es = 2/0,E + (1 + A y ,s) 



-x Es sin cf) ViE + y Es cos <j> ViE ] 



(A5) 

Note that our transformation uses 6 parameters (as does that of lAK04l ). not just the normal four parameters of reference 
posit ion, position angle, and scale, but we can also set <p y E — X y E — to have a four-parameter transformation (see 
Add. [TO. 

In the fit, we use Eqs lA4l and lA5l to calculate predicted positions, and fit simultaneously for E E and S s by minimizing, 



* 2 = E E 

s=l S=l 



( %E 

V a x,E 



%Es 



^Ve_ 



VEs 



0~y.Es 



(A6) 



This expression is very similar to the one for the combination of the exposures (Eq. IA2|) , the only difference being that 
for the combination the positional parameters (the average positions) were independent of exposure number, while 
here there is a dependence on the time of each epoch (Eq. IA4p . Because of this similarity, our implementation uses 
the same fitting routine for both steps. 

Note that the above method departs from that used by lKvKA02t which solved first for £ E , then S s , and iterated. 
We found that while our iterative solution was not biased, direct minimization was more precise, could more correctly 
disentangle the effects of stellar proper motions and uncertainties in £ E , and found more easily the global y 2 minimum. 

Our fit automatically yields a global x 2 > but, like for the exposure combination, we can also calculate \ 2 values for 
each individual star (Table [4]) and each individual epoch (Table [5]), and these are useful in identifying problems. We 
stress again, however, that these are not rigorous \ 2 values, because the exposure parameters £ E are common to all of 
the stars, and the \ 2 f° r a single star does not take this into account (similarly, S s are common to all of the epochs, 
and the \ 2 f° r a single epoch does not take this into account). As for the combinations, this affects the brightest stars 
in particular; from Table [4j one indeed sees that stars A, B, 118, and 119 have x 2 e d substantially below unity. 

B. ALTERNATE SCHEMES FOR THE PARALLAX MEASUREMENT 

We made a number of choices in determining the positions of stars in individual exposures and in combining these 
positions measurements. While we had reasons to prefer some schemes over others, we wished to verify that our final 
measurements did not depend on the detailed choices that we made. Therefore, we present the consequences of different 
ways of measuring positions (§ IB.ll) . estimating uncertainties (? IB.2[) . identifying "bad" measurements IB.3|) . and 
combining exposures into averages and fitting these to an astrometric model (§ IB.4I) . We also discuss a few statistical 
tests of the robustness of our results f? IB.5|) and end with a summary of the possible sources of systematic error IB.6|) . 

We note that all of the steps in the analysis (§ 12.21 § [4] and App. [A]) were performed independently by two of us 
(DLK and MHvK), using separate routines. We cross-checked the results at many points in the analysis, ensuring 
results were identical. 



B.l. Alternate Position Measurements 

At the start of the project, we wished to see if the measurement scheme of AK04 (Appendix C) was optimal. We did 
a number of experiments, and settled on a slightly different scheme, described in § 12.21 Overall, we found that, for our 
project, the choice of scheme had little influence: differences in t he fina l parallax were less than 0.15 mas, and we find 
the same parallax also if we simply use the routines provided bv lAK04l and do careful rejection of outliers (App. IB.3|) . 
However, the differences may be important for brighter sources, and hence we briefly descri be our experiments below. 

First, we varied the way pixels are weighted, taking into account only Poisson noise, as in I AK041 . also read-noise, or 
even a small, few percent "flat-field" error (proportional to the flux). We found that using just read-out and Poisson 
noise provided superior \ 2 f° r t ne combination of exposures. We also tried determining weights based on the ePSF 
model instead of the observed counts. Here, using observed counts has the disadvantage that statistical fluctuation 
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towards low (high) counts get too high (low) a weight), while using a model has the problem that, due to focus changes, 
etc., the model may not be precise. The results, however, where not significantly differe nt. 

Second, we experimented with the fit region, trying the 5x5 pixel box used by I AKO j . smaller and larger boxes, as 
well as circular regions with a range of radii. Overall, we found that circular regions gave slightly smaller residuals 
in combining exposures. We picked a radius of 2.56 pixels to minimize the dependence of the number of pixels in the 
region on centroid position (the average number included is 20.6, with a standa rd deviation of 0.7). 

Third, we tried determining the sky from annuli around sources, like IAK04I do, including it in the fitting process, 
and fixing it globally to the header value mdrizsky. We found only slight differences; our choice of a constant sky 
level was based on a slightly smaller x 2 ■ 

Fourth, we tried fitting for the ePSF amplitude as well as fixing these amplitudes using the average photometry and 
an aperture correction for each exposure. It had minimal effect: using fixed amplitudes gave % 2 worse by a few percent 
in the combination and astrometry. 



B.2. Alternate Uncertainty Schemes 

We determined uncertainties in the positions from the \ 2 fitting. To verify these, we simulated the detection of 10,000 
stars of a range of brightnesses, and calculated standard deviations. This simulation included the effects of photon 
noise and read noise (which dominate for the fainter stars), but did not include any mismatches in the ePSF (which 
will dominate for brighter stars). We find that the la uncertainty in each coordinate can be described approximately 

by 

1/2 



cr(rai) 



pixel, (Bl) 



where m\ = —2.5 logiV with N the counts in a 5 x 5 pixel box (or, equivalently, the ePSF amplitude from the fit), and 
C\ ~ 1.14 and C'2 — 36 are constants determined from the simulation. This relation roughly matches expectations, 
since the uncertainty should scale as the full-width at half-maximum (FWHM) divided by the signal-to-noise ratio 
(SNR), where the SNR has contributions from Poisson noise (ex y/~N; first term, which dominates for bright stars) and 
read-noise (cx 7V/RN; second term, which dominates for faint stars; here, RN is the read-noise in the pixels covered by 
the PSF). With a gain of unity for flat- fielded images and a FWHM close to 1 pixel, one expects C% ~ 1, as we find. 
Similarly, from the read- noise of ~6 DN per pixel, and considering that 20 pixels are used, one roughly reproduces Ci. 
In Fig. [5J we compare the predictions with the standard deviations found after combination of the exposures (Eg. IA3[) . 
and find good agreement. 

For the brightest stars, one might expect that systematic effects to star t to domi nate, but from neither simulations 
nor observations do w e find evidence for this (Fig.[5J), in contrast to what lKvKA02l found for WFPC2 data, and what 
is found for HRC by IAK04 For the observations, however, one has to keep in mind that, for the brightest stars, 
the standard deviations underestimate the measurement uncertainty, since part of the measurement error has been 
absorbed in the transformations required to put the exposures on a single reference frame (App. |Aj). Hence, any 
systematic error might be hidden in the transformations. 

To test the effect of possible systematic errors, we added, in quadrature, constant terms of up to 0.01 pixel to 
our input uncertainties, thus reducing the weight carried by the brightest stars. As expected, the overall x 2 f° r the 
exposure combination and the astrometry decreases, from Xred.comb = x\l J2e ^dof.B = 1-21 and X^ed.par = 1-0? 
to 0.89 and 0.81, respectively. The change in parallax, however, is only 0.005 mas. We can use these results to estimate 
the size of the systematic errors, by determining for w hat add itional systematic error one obtains X 2 e d — We- find 
this is at about 0.007 pix, consistent with the results of lAK04 



B.3. Alternate Rejection Schemes 

During our analysis, we found that a very important part in obtaining a good registration is the rejection of 
measurements that are biased by a bad pixel or a cosmic ray. The scheme we used was meant to be as objective 
as possible, being based on identifying bad pixels in an entirely automatic way. Yet, we still needed to reject some 
additional measurements (§ I2.2j) . As an alternative, we identified bad measurements from large deviations of the 
positions and/or instrumental magnitudes of single observations from the average. We did not employ a formal 
rejection criterion but examined all of the data manually, iteratively registering the exposures and removing outliers. 
For the brightest four stars, very few exposures were rejected, while for the fainter ones we included on average 
7 exposures per epoch. With this set, we derived a parallax very similar to our best value. 



B.4. Alternate Combination Schemes 

We combine our individual measurements into averages using a six-parameter transformation between exposures, 
and again use a si x-parameter transformation in our astrometric solution. Although these choices are grounded in 
the work of [AK04, we tried some alternatives: (i) combining exposures with two or four-parameter transformations; 
(ii) tying the epochs using a four-parameter transformation; and (iii) instead of merging the exposures, treating each 
exposure as a separate epoch, and fit for 64 transformations £s simultaneously. 

We found that a different choice of combining the exposures did not change the final parallax significantly (i.e., by 
more than 0.1 mas), although it did affect the quality of the fit: with fewer parameters, the x 2 values were significantly 
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higher. Using only 4 parameters for the astrometric tie has a more drastic effect on the quality of the fit. The parallax, 
however, decreased by only 0.1 mas compared to our main result. 

When we treated all exposures as separate epochs, we found a parallax identical to within 0.01 mas with that 
derived using our two-stage approach. The fit was formally less good than our two-stage result, which may reflect an 
underestimate of the uncertainties of bright stars and/or remaining position measurements contaminated by cosmic 
rays. We made a number of experiments to address these issues, but, fortunately, the effects of all of them on the 
parallax is small, < 0.1 mas. 

B.5. Statistical Tests 

As a statistical measure of the robustness of our result, we tried to "jackknife" our data, removing either one star 
or one epoch from the fit. Removing a single star does not change the fit appreciably, showing that no single star 
dominates the fit. On average, there was no change in the parallax, and the dispersion was only 0.03 mas. The 
dispersion is dominated by one star, B: with B omitted, the parallax decreased by 0.1 mas (note that B is also the 
most discrepant point in Fig. |4])- Given how small this change is compared to the statistical uncertainty, we are 
satisfied with this test. 

Removing an epoch had a more dramatic effect. As can be seen from Fig.[6l some epochs fit significantly better than 
others. Specifically, removing epoch 5 improved the quality of the fit (Xred par = 0-93 for the overall transformation), 
and increased the parallax to 3.32 ± 0.96 mas. In contrast, removing epoch 8 gave a less significant change in the 
quality of the fit (x;? cd par = 1.03), but a larger change in parallax, to 1.67 ± 1.01 mas. These were the most extreme 
examples: the mean shift was 0.08 mas, and the standard deviation was 0.5 mas (0.3 mas excluding the epoch 8 result). 
Part of the problem with removing epoch 8 seems to be a covariance with the proper motion: while removing the other 
epochs changed the proper motion by < la, removing epoch 8 changed [i a to —91.9 ± 1.5 mas yr _1 (the correlation 
coefficients between 7r and fi a or fis are significant; see Table [6]). 

Finally, we also removed individual exposures when we fit using all of the exposures without combination (App. lB~4|) . 
Again, the mean of the parallaxes was identical to that in Table El and the standar d deviation was 0.1 mas. This 
leads to an uncertainty of 0.9 mas for the parallax ( following lEfron fc Tibshiranilll993t ), which is the same as what we 
derived from error propagation. 

B.6. Estimate of the Systematic Uncertainties 

The biggest source of systematic uncertainties that could affect our results are the epochs that we include. While 
the other choices, as des cribed above, led to parallax differences of < 0.1 mas, the jackknife tests on epochs gave 
larger changes (App. IB.5[) . Overall, we believe our result from Table [5] to be reliable, but we should add a systematic 
uncertainty to the statistical uncertainty quoted there. The magnitude of this uncertainty should be similar to the 
variation in the parallax that we see, or 0.4 mas. This may be an overestimate, as our exposure jackknife tests indicated 
that the bias was small, but we prefer to err on the conservative side. A small number of measurements in the future 
that will have a long enough time baseline to make the proper motion better determined should reduce this uncertainty 
significantly. 
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